Frequency dependent polarizability of small metallic grains 
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We study the dynamic electronic polarizability of a single nano-scale spherical metallic grain 
using quantum mechanical approach. We introduce the model for interacting electrons bound in the 
grain allowing us numerically to calculate the frequency dependence of the polarizability of grains 
of different sizes. We show that within this model the main resonance peak corresponding to the 
surface plasmon mode is blue-shifted and some minor secondary resonances above and below the 
main peak exist. We study the behavior of blue shift as a function of grain size and compare our 
findings with the classical polarizability and with other results in the literature. 

PACS numbers: 73.20.Mf, 73.22.Lp, 78.67.Bf 

INTRODUCTION 

Great efforts in contemporary materials science research focus on properties of granular materials [lj ■ The interest is 
motivated by the fact that granular arrays can be treated as artificial solids with programmable electronic properties. 
The ease of adjusting electronic properties of granular materials is one of their most attractive assets for fundamental 
studies of disordered solids and for targeted applications in nanotechnology. The parameters of granular materials in 
many ways are determined by the properties of their building blocks - grains. 

In this paper we study the dynamic electronic polarizability of a single nano-scale spherical metallic grain. This 
research is important because on one hand confinement effects can drastically modify the dynamic polarizability 
of metallic grains and on the other hand the fundamental question related to the description of crossover between 
microscopic and macroscopic behavior of the grains exists. 

The electronic polarizability was first quantum mechanically studied in Ref. 0]. However, the influence of screening 
effects on static polarizability was first investigated in Ref. [3j . Later the dynamic polarizability of very small metallic 
particles was considered within the jellium-background model ■ This model was also used to study the electronic 
structure and polarizability of metallic nanoshells (s | . The electronic polarizability was also studied in applications 
for optical properties of small metallic grains More recently the nonperturbative supersymmetry technique 

was used to consider the effects of disorder on properties of electronic polarizability [3] ■ 

Here we study the dynamic electronic polarizability of a single metallic grain using quantum mechanical approach. 
We introduce the model for interacting electrons bound in the grain allowing us numerically to calculate the frequency 
dependence of the polarizability of grains of different sizes. 

MAIN RESULTS 

We numerically calculate the dynamic polarizability a(uj) of grains of different sizes. The experimentally measurable 
quantity is the photoabsorption cross section a(u>). It is related to the imaginary part of the polarizability a(oj) as 
follows 

a(ui) = Imev(cj), (1) 

c 

where c is the speed of light. The photoabsorption cross section is the main result of this work. It is plotted in 
Figs. [TJ and [2] for two grains of different sizes (the number of electrons N = 16304 and N = 1956, corresponding to 
the grain sizes of 5.28 nm and 2.6 nm, respectively). 

Our calculations were done at zero temperature. Therefore, due to Fermi-Dirac statistics all electron states are 
occupied below the Fermi level and empty above it. As a result, the total number of electrons N, that determines 
the size of the grain, is related to the Fermi energy level. This number can be found by calculating the total number 
of electron states below the Fermi level. The Fermi energy can be parameterized by dimensionless quantity aj? as 
Ep = (h 2 /2mR 2 )a F , where R is a grain radius. In order to estimate the grain size we need to know the total number 
of electrons N and the electron density r s . In terms of these parameters the grain size is given by R/as — A 1 / 3 ^, 
where as is Bohr radius. All grain size estimates were done for sodium (r s = 3.93). 

The behavior of dynamic polarizability a(ui) presented in Figs. Q] and [5] revels several features: i) the main resonance 
peak is slightly shifted (with respect to the classical result) to higher frequency; ii) some smaller secondary resonances 
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FIG. 1: (Color online) Imaginary part of dynamic electronic polarizability a(ui) of the metallic grain (in units of R 3 , R being 
the radius of a grain) vs dimensionless frequency uj/lj c i with lj c i — cj p /\/3 being the classical resonance frequency, where uj p is 
the plasma resonance frequency. The blue solid line represents the numerical solution for polarizability a(ui) corresponding to 
the number of electrons N = 16304, R — 5.28 nm(QF = 50, see Eq. (|22[) ). The red dashed line represents the classical Drude 
polarizability a c i{io), Eq. ([3]), for comparison. 

above and below the main resonance frequency exist. Both of these features have a tendency to decrease as the grain 
size increases which is consistent with the fact that macroscopic grain has a fully classical behavior. 

The main resonance peak in Figs. [Hand [2] corresponds to the surface plasmon excitation. The blue-shift of the main 
resonance peak in comparison with the classical Mie resonance (the red dash line in Figs. [T] and [2]) can be understood 
as follows: In a quantum mechanical problem due to Fermi-Dirac statistics electron transitions may occur only from an 
occupied state to an empty state. For higher external frequency the total number of unrestricted electron transitions 
becomes higher. Therefore in a quantum mechanical problem if the external frequency equals to the classical plasmon 
frequency, u> c i = oj p /\/3 in Figs. Q] and [2] the resonance amplitude is suppressed because of a lack of possible electron 
transitions. Hence the main resonance peak has a tendency to move to a higher frequency. 

The surface plasmon excitations were also studied quantum mechanically in Ref. [15| , and using extended classical 
Mie theory in Ref. 10], where it was shown that the frequency of a surface plasmon for very small grains is blue- 
shifted. This is in agreement with our result. Several experimental works have also observed the blue-shifted surface 
plasmon resonance peak [la . Il7l ] . 

A different result was obtained in Refs. [Hill: where the main resonance peak was red-shifted. This discrepancy 
may be explained by the difference between the models used. In particular, in Ref. [4| the electronic polarizability 
was studied within the spherical jellium-background model. The red-shift in this model is explained by the electron 
spill-out effect. Here we use "electron in a box" approximation. In this model the spill-out effect is prevented by the 
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FIG. 2: (Color online) Imaginary part of dynamic electronic polarizability a(ui) of the metallic grain (in the units of R 3 , R 
being the radius of a grain) vs dimensionless frequency uj/ujci with uj c i = oj p /\/3 being the classical resonance frequency. The 
blue solid line represents the numerically solution for polarizability a(cj) corresponding to the number of electrons N = 1956, 
R = 2.6 nm (oif = 25, see Eq. (|22|) ). The red dashed line represents the classical Drude polarizability a c i(ui), Eq. Q, for 
comparison. 



infinitely deep potential well. 

The width of the plasmon resonance in Figs. [1] and [2] is determined by the artificially introduced imaginary part of 
the frequency u> to avoid infinite values of polarizability cv(w) at the resonances. This width does not correspond to a 
real width of the photoabsorbtion cross section <t(ui) in Eq. ([!}. We discuss this issue in more details in the "solution" 
section below. 

The secondary resonances shown in Figs. [T] and [5] are significantly lower than the main peak. These resonances can 
be divided into two groups: i) low frequency single particle excitations and ii) high frequency resonances corresponding 
to the transitions between the states with higher orbital quantum numbers. Both of these types of excitations become 
less pronounced with increasing grain size. We notice that at low energies the transitions may occur only in the 
immediate vicinity to the Fermi surface. Therefore the low frequency peaks in Figs. Q] and [2] correspond to the 
resonant transitions between the states which arc located just below and just above the Fermi surface. 

We now comment on the influence of disorder in the grain on the resonance peaks. First, we notice, that our 
consideration is based on the fact, that the disorder is weak. Thus, the conducting electrons are dclocalized in 
our model. The characteristic energy scale for disorder is the inverse scattering time 1/r. In the limit of weak 
disorder the electron propagation within the grain is almost ballistic and therefore 1/r is of the order of Thouless 
energy, 1/r ~ Exh- However, both of these energy scales are smaller than the plasma frequency, uj p , which is the 
characteristic scale for the resonances. Therefore, we expect that for energies less than lo p the disorder is irrelevant 
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FIG. 3: (Color online) Blue shift of the surface plasmon frequency (in the units of the classic surface plasmon frequency 
LJd = Lo p /y/3, with u) p being the plasma frequency) as a function of the dimensionless grain radius, R/Xf, where R is the grain 
radius and Xf is a Fermi wavelength. Red dots represent the numerical calculations. The blue solid line represents the fitting 
of the numerical calculations with the function (R/Xp)' 1 



and can not modify the picture qualitatively. 

However, the presence of disorder may lead to the broadening and suppression of all of the resonance peaks. The 
broadening effects are beyond our present consideration. The disorder caused suppression might make the secondary 
resonances hard to detect, as they are small compared to the main resonance peak. 

We also mention that our consideration can be generalized for the description of slightly non-spherical grains. In 
this case the non ideal grain shape can be studied within the perturbation theory. It will result in the changes of the 
energy levels structure. In particular, the degenerate energy levels will be split. That will slightly affect the positions 
of the resonance peaks. 

To conclude the analysis of the polarizabilities we compare our findings for a(u>) with the known result for classical 



polarizability a c /(w), 21 [ 



e(w) + 2 

where R is the radius of a grain and e(w) is the grain dielectric constant; within the Drude model e(uj) = 1 — ^j^y, 
where w p = ((4/ 3)nh') 1 ^ 2 ev p is the plasma resonance frequency and S is the damping factor. The imaginary part of 
the polarizability a c i (w) in Eq. ^ is given by 

lma c i(oj') = R 3 - — 2~^2' ( 3 ) 

where we introduce the dimensionless parameter 7 = 5/(lu p /\/3) and the dimensionless frequency oj' = uj/ (oj p /\/3). 
The right hand side of Eq. ([3]) has a single resonance peak at frequency oj' = 1 (or u> = cj p /\/3) corresponding to the 
surface plasmon mode excitation. We plotted Ima C !(u) in Figs. Q] and [2] for comparison with our results. 

The last important question to discuss is the dependence of the surface plasmon frequency on the grain size. This 
dependence is shown in Fig. [3j The calculated red dots were then fitted with the function that scales as (R/Xf)^ 1 - 
On can clearly see that the blue shift has the tendency to decrease as the grain size increases. 

Now we turn to the description of our model and the derivation of electronic polarizability a(u) in Eq. JT]). 

THE MODEL 

A single spherical metallic grain is described by the following Hamiltonian 

H = H + H int . (4a) 
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Here 

H 



dW(r) f-J- + U(t)] ip(r), (4b) 
\ 2m ) 



is the Hamiltonian of the non-interacting electrons trapped in the infinitely deep spherical potential well U(r): 

This form of the potential allows us simple analytical expansion in eigenfunctions of the Hamiltonian Hq, which in its 
turn allows us to use the many-body approach. The boundary conditions in Eq.0c]are valid if the number of electrons 
near the grain surface is small compared to the total number of electrons in the grain, i. e. R/Xf 3> 1, where R is 
the grain radius and Xf is Fermi wavelength. In Eq. (|4bl) ?/^( r ) [VK r )] is the creation [annihilation] operator and the 
integral is performed over the grain volume. 

The second term in the r. h. s. of Eq. (|4a|) is responsible for Coulomb interaction 



H int = J d^^r^^r')——^)^'). (4d) 

The common method of studying the electronic properties of such systems is based on a mean-field approach. 
This approach reduces the many-body problem to a single particle problem with renormalized effective potential, 
where the effective potential contains all of the corrections due to the interaction between the electrons. However, in 
the present manuscript we use a different approach. Using the many-body Hamiltonian, Eq. (|4al) . we consider the 
Coulomb interaction, Eq. f|4d[) . perturbatively. This allows us to tackle the problem without reducing it to a single 
particle problem with renormalized potential U(r). We perform our calculations in the basis of the eigenfunctions 
of the non- interacting Hamiltonian Hq, Eq. (|4b|) and take into account interaction effects, Eq. (|4dp . summing the 
appropriate series of Feynman diagrams. 

Our consideration is valid for metallic grains meaning that the grain itself is a good conductor. This means that 
the grain size R is much larger than the Fermi wavelength A^. In addition, the disorder within the grain is weak such 
that electrons within the grain are delocalized. If these conditions are valid there is a small parameter in the problem 
- the invers grain conductance 1/g, where g is the grain conductance. The perturbation theory is built then on the 
expansion with respect to this small parameter. 

To calculate the dynamic polarizability a(w) we study the electronic response of a single metallic grain to the 
external electric potential $ ea: *(r)e _4Ut . The polarizability is defined as a coefficient between the external electric field 
E and the induced dipole moment, p(w) = a(w)E. Thus, the problem is reduced to find the electric dipole moment 
p(aj) excited by the external field 

p(w) = / d 3 rrp m (w,r), (5) 

where p m (u>,r) is the field-induced electronic charge density. Within the linear response theory p m (w,r) is given by 
the following expression, see e. g. Ref. [3] 

p in (u, r) = - J dV$(w, rOx(w, r, r'). (6) 

Here x( w i r i r ') is the electron density-density correlation function and $(w,r) is the total electric potential, which 
is a sum of the external potential $ ext (w,r) and the induced potential $ m (w,r), $(w,r) = $> ext (uj,r) + $ m (w,r). 
Equation ([6]) corresponds to the random phase approximation [19j . 

In the quasi-static limit the potential $(w,r) can be found by solving Poisson equation V 2 $(u;,r) = — 47rp m (w,r). 
Using this expression in Eq. ([6]) we obtain the following integral equation for the potential 

V 2 $(w,r)=47r J dV$(w,r')x(w,r,r'). (7) 

In Eq. ([7]) the electron density-density correlation function x( w : r : r ') is expressed in terms of electron Green's functions 
G(e,r',r) as x(u;,r,r') = — TY"], G(e + u, r, r')G(e, r' , r). Using the properties of Green's functions x(w,r, r') can be 
written in the form, see e. g. Ref. [20| 

X (^,r,r') = 2 e ^/Ml^l C(r )^(/ )ft( /) i , (r)i (8 ) 



6 



where f(e n ) is the Fermi-Dirac distribution function with e„ being the energy level of the quantum system and ip n (r) 
is the eigenfunction of the non-interacting Hamiltonian with n being the full set of quantum numbers that describe 
the state. Below we consider the solution of Eq. (J7]) in details. 



SOLUTION 



To solve Eq. ([7} we introduce the eigenfunctions ip n im and corresponding eigenvalues ei n for electrons in the spherical 
potential well U(r) of radius R 

Tpnim = Pi n ji{ot ln r/R)Yi m (6,(t>), ti n = af n /2mR 2 , (9) 

where a/„ is the n-th zero of the spherical Bessel function ji(x), Yi m (9,(f>) is the spherical harmonic, and /3/„ = 
(2 1 / 2 /i? 3 / 2 )j i '(a; n ) is the normalization coefficient with j[{ai n ) being the derivative of the Bessel function taken at 

The potential $(w,r) in Eq. (0 consists of two parts: applied external potential and induced dipole part 0]. 
Choosing the direction of the z axis along the external electric field E u we obtain 

$(w,r) = - E u r cos 9 + <5> m {uj,r). (10) 

Here the induced electric field $ m (w,r) vanishes at infinity, \im\ T \^ f00 'S/^ m {uj,r) = 0. In addition, the potential 
$'™(u,r) and its derivative are continuous everywhere including the grain surface. Outside the grain the charge 
density p m (w,r) is zero therefore the potential § m (ui,r) must satisfy the Laplace's equation 

V 2 $ m (w,r) = 0, |r| > R. (11) 

Thus, the potential outside the grain can be found as a classical dipole potential 

**{ U ,T)= C °( u)E 'f C ° S9 , (12) 

where cq(u)) is yet unknown function of frequency u>. 
Inside the grain one can seek the solution in the form 

$ m (w,r) = E u [c (u>)r + Rf(u>,r)]cose, (13) 

where f(tj,r) is some unknown function of frequency u) and distance r. Introducing the function g(u,r) as 

p^(uj,r) = ^g(u,r) C o S 9, (14) 

47T 

we can rewrite Eq. ([6]) as follows 

,R 

g(u,r)= / dss 2 q(uj,r,s)[(c (u;) - l)s/R+ f(u>,s)}, (15) 



with the kernel q(u>,r, s) being 

q(u, r, s) 
Ait 

Since we are interested in the dipole response of the grain we expand the functions /(w, r) and g(u>, r) in Eq. ([15 
using the orthogonal set (f> n (r) — j\{ai n r / R) 



= / cos6dcj)d6 / cos0'dcj>'d9'Y w (9)Y w (9')x(uJ,r,s). (16) 



OO 

n— 1 

g(u,r) = E^ItFMt)" (17) 



n=l 
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Thus, the problem of calculating the induced potential $ m (w, r) inside the grain is reduced to the problem of calcu- 
lating the frequency dependent coefficients c n (uj) in Eq. (|17p . Using these coefficients the dynamic polarizability a(ui) 
of a grain can be written as follows 

a(u) = c q (uj)R 3 , (18) 

with R 3 = a c i(0) being the static classical polarizability. The coefficient cq(uj) can be considered as a dimensionless 
dynamic polarizability and can be found using the boundary conditions 

co(u>) = --V c n (uj)j[(a ln )a ln . (19) 

6 n=l 

To calculate the frequency dependent coefficients c n (uS) in Eq. (fT9| we substitute Eqs. ([T7|) into Eq. (p~5)) to obtain 
the system of linear equations for coefficients c n (u>) 



where 



A mn (ui) c n (uj) = a m (w), (20a) 



/ ,-, 3 X F X F / fim — /(?-i)n 2 i D i /j \riin ^ 

am ( w ) = TwfT V 5" ~ Z B m('. «i, n 2 )C (Z, rii, ra 2 ) 

2V27T V a B R l ^ W ~ + e (i-l)n 2 7ini7(i-l)n 2 

Jinj ^ /(Z+l)n 2 i + 1 D 2 ^2/i \\ / nnk \ 

-B m {l,ni,n 2 )C (t,ni,n 2 ) , (20b) 



w ' - £lm + e (f+l)n 2 7/ni7(i+l)n 2 



. , m 3 Xf Xf \^ / //ni f(l-l)n 2 I D l /, i D l/, s. 

A mn {u ) = p ~ T — B m {l ) n u n2)B n {l ) ni,n 2 ) 

2V2n \ a B R W - e lni + e ( /_i)„ 2 7toi7(*-i)n 2 



flm — f(l+l)n 2 Z + l 2 2 



w' - £;m + e (;+i)n 2 7(m7(/+i)n 2 



B m (l,ni,n 2 )B n (l,ni,n2) 



-^a m {^')ainj'i(o: ln ) - 5 mn a\ n ^i n . (20c) 

Here as = l/rne 2 is the Bohr radius and we introduce the notation -fi n = ^ [j'i{&in)] 2 , and the dimensionless frequency 
uj' = uj/(u} p / \/S), with ujp = (4/3)TTi i e 2 v F being the plasma frequency. 

The quantities m, n 2 ), B^l, rii, n 2 ), C 1 (/,ni,n 2 ), and C 2 (/,ni,n 2 ) in Eqs. (|20b[) and (|20c[) are defined as 

s m 2) (') n i> n 2) = / rfra 2 ji(ai„a;)j;(a; ni a;)j; T i(a( iTl ) n2 x), 
Jo 

C 1 ^^/,^!,^) = / da;x 3 j;(a;„ 1 x)7 ;Tl (Q!( ZTl) „ 2 a;), (21) 
Jo 

where the upper index 1 (2) corresponds to — (+). The Fermi-Dirac functions /;„ in Eqs. (|20b[) and (|20c|) are taken 
at zero temperature and defined as 

fm = {]: ain t aF ( 22 ) 

[ 0, ai n > OF, 

where the dimensionless parameter ckf defines the Fermi energy E F = {h 2 /2mR 2 )a F . Using Eq. ([22]) for the distri- 
bution function //„ one can calculate the number of the electrons in the conducting band as follows 



AT = 2 z ^(2Z + l)/ ta . (23) 



For numerical calculations we used the complex frequency uj' + v~f with 7' = 0.006 being the dimensionless damping 
factor. The damping factor 7' was introduced to avoid infinite values of the polarizability at the resonances and does 
not correspond to the real plasmon resonance broadening. The analysis of the resonance broadening is beyond the 
scope of the present work. 

Using Eqs. (fT5)) - (|20cp we numerically calculate the dynamic polarizability a(uj) in Eq. (fTS)) for grains of different 
sizes. The final results are shown in Fig. [1] and Fig. [5] 
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CONCLUSIONS 



We studied the dynamic polarizability of spherical metallic grains using quantum mechanical treatment. We 
numerically investigated the frequency behavior of polarizability for relatively large grains. We showed that the main 
resonance peak corresponding to the surface plasmon mode is blue-shifted and some minor secondary resonances above 
and below the main peak exist. We studied the dependence of blue shift as a function of grain size and compared our 
results with the classical polarizability. 

We thank G. Weick for helpful discussions. This research was supported by an award from Research Corporation 
for Science Advancement. 
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